CMB Maximum Temperature Asymmetry Axis: Alignment with Other Cosmic 

Asymmetries 
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We use a global pixel based estimator to identify the axis of the residual Maximum Temperature 
Asymmetry (MTA) (after the dipole subtraction) of the WMAP 7 year Internal Linear Combination 
(ILC) CMB temperature sky map. The estimator is based on considering the temperature differences 
between opposite pixels in the sky at various angular resolutions (4° — 15° and selecting the axis 
that maximizes this difference. We consider three large scale Healpix resolutions (N 3 id e = 16 (3.7°), 
N S id e = 8 (7.3°) and N s ide = 4 (14.7°)). We compare the direction and magnitude of this asymmetry 
with three other cosmic asymmetry axes (a dipole, Dark Energy Dipole and Dark Flow) and find 
that the four asymmetry axes are abnormally close to each other. We compare the observed MTA 
axis with the corresponding MTA axes of 10 4 gaussian isotropic simulated ILC maps (based on 
ACDM). The fraction of simulated ILC maps that reproduces the observed magnitude of the MTA 
asymmetry and alignment with the observed a dipole is in the range of 0.1% — 0.5% (depending on 
the resolution chosen for the CMB map). The corresponding magnitude+alignment probabilities 
with the other two asymmetry axes (Dark Energy Dipole and Dark Flow) are at the level of about 
1%. We propose Extended Topological Quintessence as a physical model qualitatively consistent 
with this coincidence of directions. 

PACS numbers: 98.80.Es,98.65.Dx,98.62.Sb 



1. INTRODUCTION 

There is observational evidence coming mainly from 
the isotropy of the Cosmic Microwave Background 
(CMB) that the Universe is isotropic on Hubble scales. 
Any anisotropy on these scales is bound to be smaller 
than about 1 part in 10 3 . This constraint combined with 
the Copernican principle (supported by kSZ data [l[) 
leads to strong support of the cosmological principle: the 
Universe is homogeneous and isotropic on Hubble scales. 

The violation of the cosmological principle is expected 
to occur at a small level even on Hubble scales due to 
small statistical fluctuations of the cosmic energy den- 
sity (matter, radiation, dark energy). Precise cosmolog- 
ical observations are in principle able to detect this de- 
viation from the cosmological principle on Hubble scales 
and compare the expected magnitude with the one an- 
ticipated based on the standard isotropic cosmological 
model. 

The lowest order deviation from isotropy, which is also 
easiest to detect, is the anisotropy that distinguishes a 
preferred cosmological axis. Such an axis is usually rea- 
sonably described by a dipole deviation from isotropy. 
An exception consists of the CMB temperature pertur- 
bations where the dipole term is dominated by our mo- 
tion with respect to the CMB frame and therefore it has 
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been removed completely from the CMB maps. This re- 
moval has also swept away any subdominant cosmological 
contribution to the dipole. However any axial cosmologi- 
cal anisotropy that is not perfectly described by a dipole 
could have left a trace after the removal of the dipole. 
The detection of this trace may be possible by using spe- 
cially designed statistical tests. 

Early hints for deviations from isotropy on Hubble 
scales have been accumulating during the last decade. 
Some of these hints may be summarized as follows (HQ: 

1. Large Scale Velocity Flows (Dark Flow): 

There are recent indications that there is a large 
scale peculiar velocity flow with amplitude larger 
than 400km/sec on scales up to 100ft. _1 Mpc (z < 
0.03) [4] with direction I ~ 282° ±11°, b ~ 
6° ± 6°. Other independent studies have also 
found large bulk velocity flows with similar direc- 
tion [5| on scales of about 100/i _1 Mpc or larger 
than 300/i _1 Mpc Q. This large scale peculiar ve- 
locity flow is known as Dark Flow. The stan- 
dard homogeneous- isotropic cosmology (ACDM), 
predicts significantly smaller amplitude and scale 
of flows than what these observations indicate. The 
deviation of these observations from ACDM predic- 
tions is more than 3er. Other studies [7] using Type 
la supernovae find a flow of a somewhat smaller 
magnitude, consistent with ACDM. Even in these 
studies however the direction of the flow is similar 
to the direction found by Refs. [4j and [6|. Thus, 
even though there is some controversy about the 
magnitude of the Dark Flow it appears that its di- 
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rection is more robust even though its directional 
la error region is probably larger than the one in- 
dicated in Ref. H (for a more conservative la error 
region see Ref. [7] ) . A possible connection of large 
scale velocity flows and cosmic acceleration may be 
found in Ref. j|. 

2. Fine Structure Constant Dipole: Quasar ab- 
sorber spectra obtained using UVES (the Ultravi- 
olet and Visual Echelle Spectrograph) on the VLT 
(Very Large Telescope) in Chile and also previous 
observations at the Keck Observatory in Hawaii 
Q have indicated that the value of the fine struc- 
ture constant at high redshifts (z s [0.2,4.2]) is 
not distributed isotropically. Rs anisotropy is well 
described by a dipole with axis directed towards 
I ~ 320° ± 11°, b ~ -11° ± 7°. The deviation of 
these observations from isotropy is 4.1a (9l. Hoj. 

3. Dark energy Dipole: A recent fit of the Type 
la distance moduli residuals (from the best fit 
ACDM) to a dipole anisotropic distribution has in- 
dicated [l(| that the angular distribution of these 
residuals is well described by a dipole analogous to 
the Fine Structure Constant dipole. Rs axis is to- 
wards I ~ 309° ± 18°, b ~ -15° ± 11° and deviates 
by only about 11° from the fine structure constant 
dipole. The deviation of these observations from 
isotropy is at the 2a level. 

Each one of the above observed deviations from isotropy 
is between 2a and 4ct. The angular proximity of the cor- 
responding anisotropy axes makes their combination even 
more unlikely in an isotropic universe where there is no 
correlation between them. In Ref. [l(| it was shown that 
the combined magnitude-l- alignment of the fine structure 
constant a and dark energy dipoles has a probability less 
than one part in 10 6 to occur in an isotropic universe 
where the two dipoles are uncorrelated. 

A physical model was proposed in Rcfs. [HI El that 
has the potential to explain the existence and the align- 
ment of the above three dipoles. The model is based 
on the existence of a topological defect (e.g. a global 
monopole) with Hubble scale core formed during a recent 
phase transition by an 0(3) symmetric scalar field non- 
minimally coupled to electromagnetism. An off-center ob- 
server with respect to the monopole center would observe 
faster accelerating expansion towards the core where the 
vacuum energy density is larger and also varying a along 
the same direction due to the variation of the scalar field 
magnitude. This model is a generalization of 'Topological 
Inflation' [121 ] and has been called Extended Topological 
Quintessence [ic| due to its non-minimal coupling to elec- 
tromagnetism. The model also has some similarities with 
texture models which have been considered as a physical 
origin of the observed Cold Spots on CMB maps [3l| . 
In contrast to Extended Topological Quintessence how- 
ever, texture models have not been considered as physi- 
cal origin for cosmic accelerated expansion of for spatial 
variation of a. 



Extended Topological Quintessence makes the follow- 
ing qualitative predictions (Toj : 

1. Large Scale Velocity Flows: Due to the stronger 
repulsive gravity in the defect core a large scale pe- 
culiar velocity flow is predicted along the axis that 
connects the off center observer and the monopole 
core. The direction of the flow is predicted to be 
away from the repelling core ('Great Repulser') and 
its scale is predicted to be the Hubble scale (the 
defect core scale). A reversal of the velocity flow 
direction is predicted for observations that go be- 
yond the defect core. As discussed in the following 
section, the direction of the observed Dark Flow 
is consistent with the directions of the fine struc- 
ture and dark energy dipoles in accordance with 
the above prediction. 

2. Correlation between values of a and pres- 
ence of strong magnetic fields: As discussed in 
Ref. [l(J, the scalar field magnitude is expected to 
depend weakly of the presence of local strong mag- 
netic fields. This magnitude variation is in turn 
expected to lead to local variations of a in cosmo- 
logical regions with large magnetic fields. 

3. Maximal large scale CMB variation towards 
the defect core: Due to the recent formation of 
the global defect, the ISW effect is expected to lead 
to large temperature differences between opposite 
directions in the sky along the direction towards the 
defect core. A large part of this temperature asym- 
metry would have been subtracted from CMB maps 
along with the dipole moment which is mainly due 
to our motion with respect to the CMB. However, 
smaller traces of this asymmetry could have sur- 
vived the dipole subtraction and may be detectable 
in large scale CMB maps. 

A wide range of large scale anomalies have been de- 
tected on CMB maps [HI, [13]. The anomalies in- 
clude an abnormal alignment and planarity of the oc- 
topole and quadrupole momen ts fTR 1 1 6ll . the existence of 
two large and deep cold spots |17H19j ]. the lack of large 
scale po wer 1201 -4231] . the even excess of the CMB power 
spectrum[24|. the hemispherical power asymmetry [2f| 
and quadrupolar dependence of the two point function 
(see Ref. [lj] for a detailed review). Recent evidence for 
mirror symmetry and antisymmetry (along different di- 
rections) has also been obtained [26|, (27| using the ILC 
WMAP7 CMB map[2|. Finally evidence for the exis- 
tence of statistically significant giant rings in the CMB 
sky has also been reported [2{|. Some of these anomalies 
appear to be related to a large scale bipolar asymmetry 
of the CMB even though there is no current quantitative 
study of a physical model than can give rise to all these 
anomalies simultaneously (see however [30M33| ). Due to 
the lack of such a model these anomalies are usually as- 
sumed to be a posteriori manifestations of expected large 
statistical fluctuations. 
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Having at our disposal a well denned physical model 
which makes specific predictions allows us to focus on 
specific aspects of CMB maps and search for signatures 
of our model. Thus, in what follows we focus on the pre- 
dicted large scale CMB anisotropy and search for the axis 
of maximal temperature asymmetry in the WMAP7 ILC 
map. In particular we consider three large scale Healpix 
pixelizations[34| of the WMAP7 ILC map and identify 
those pairs of opposite pixels in the sky that correspond 
to Maximum Temperature Difference. We compare the 
magnitude of this Maximum Temperature Asymmetry 
(MTA) with that expected from an isotropic model us- 
ing Gaussian simulated CMB maps. We also compare 
the direction of the MTA with the direction of the other 
observed cosmic asymmetry axes (Dark Flow, Dark En- 
ergy Dipole and a Dipole). We find the likelihood that 
the observed magnitude and alignment would occur by 
chance in an isotropic model with no correlation between 
the CMB and the other observables. 

The structure of this paper is the following: In the next 
section we describe in some detail the method for identi- 
fying the MTA magnitude and direction in the WMAP7 
ILC map. We also show the resulting magnitude and di- 
rection and also its alignment with the other observed 
axes. We then compare the observed magnitude and 
alignment with those obtained by 10 4 Gaussian simulated 
ILC maps based on ACDM. In section III we discuss the 
implications of our results and point out the next steps 
of this research program. 



2. METHOD-RESULTS 

The subtraction of the dipole moment from the CMB 
maps removes along with the dominant Doppler compo- 
nent any cosmological signal that may happen to have a 
dipole anisotropy. Such a signal is expected to emerge 
in the context of extended topological quintessence as 
discussed in the Introduction (see also Refs. [UJ, In 
addition to the dipole however, an off center observer will 
also experience axial anisotropies correspond ing to higher 
moments although at a smaller magnitude [35l [36j . De- 
pending on the dynamics and the geometry of the form- 
ing defect these higher moment asymmetries may have 
detectable magnitude. Such asymmetry could manifest 
itself as maximized temperature difference between op- 
posite pixels in the CMB sky. In order to obtain the 
direction and magnitude of such residual MTA we use 
the following steps applied on the WMAP7 foreground 
reduced ILC map pixelized according to Healpix. In or- 
der to minimize foreground contamination we focus on 
large angular scales ( N si d e = 4 (pixel size about 14.7°), 
Nside = 8 (pixel size about 7.3°), N si d e — 16 (pixel size 
about 3.7°)) 

1. Construct a Temperature Difference Map (TDM) 
obtained by assigning to each pixel a number equal 
to the difference between its temperature value and 



the value of the temperature of the opposite pixel 
in the sky. Thus we have 

D-{f H ) = {T{n i )-T{-h i )) , (2.1) 

where hi is the direction of the i th Healpix pixel. 
A similar estimator was considered in Ref. [33] in 
an effort to test statistical isotropy of CMB maps. 
In the context of the Healpix pixelization the op- 
posite pixel is always simply defined and identified. 
By construction, opposite pixels of the TDM are 
assigned to opposite values. 

2. In the TDM we select the pixel D~ aa ,(nfe) with the 
maximum absolute value. This pixel along with the 
pixel located opposite to it defines the axis of MTA. 
If the dipole had not been subtracted the MTA axis 
would be almost identical to the dipole axis. Thus 
the MTA axis is the residual asymmetry axis after 
the subtraction of the dipole. 

3. The direction of the MTA is then compared with 
the directions of other cosmic asymmetry axes (a 
dipole, Dark Energy dipole and Dark Flow) and the 
corresponding angular differences are identified. 

4. The magnitude and direction of the MTA are com- 
pared with a large number of ACDM simulated ILC 
maps [38j | and we evaluate the likelihood to obtain 
the observed MTA magnitude (or larger) in the con- 
text of ACDM. The likelihood of obtaining the an- 
gular differences (or smaller) with the other cosmic 
asymmetry axes in the context of ACDM is also 
evaluated. 

The WMAP7 ILC maps using Healpix pixelizations 
with Nside = 4, 8, 16 (corresponding to a pixel size of 
V47t/ (12iV^ de ) rad i.e. about 14.7°, 7.3° and 3.7° re- 
spectively) are shown in Figure [T] along with the MTA 
pixels. The original map has N S id e = 512. The proxim- 
ity of the MTA axis with one of the Cold Spots center is 
evident. 

In Table U we show the directions in galactic coordi- 
nates of the four cosmic asymmetry axes. In Figure [2] 
we show these directions in a Mollweide projection. The 
filled contours around each direction correspond to the 
la error regions. In Table [TTI we show the corresponding 
angular separations for each pair. 

The cumulative probability for obtaining a given value 
of MTA or larger, may be obtained usinglO 4 simulated 
statistically isotropic ACDM ILC maps 38| . The result 
is shown in Figure [3] for each one of the three angular 
resolutions (pixel sizes) considered. We used the pub- 
licly available ACDM simulated ILC maps of Ref. [Hj]. 
The observed value of the MTA magnitude is indicated 
by an arrow. The probability to obtain the observed 
magnitude of MTA (or larger) in the context of ACDM 
varies between 16% and 7% depending on the ILC map 
angular resolution. This result by itself does not indicate 
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FIG. 1: Degraded temperature maps obtained from the 7 years ILC CMB map with N a id B 
the maximum temperature differences direction. 



4, 8, 16. The white dots define 



any statistically significant deviation from ACDM predic- 
tions. Perhaps, this is the main reason that this simple 
statistic has been largely ignored by previous studies (see 
however Ref. [37j )■ However, the statistic becomes more 
interesting when the proximity of the direction of the 
MTA to other cosmic asymmetry axes is considered. 



In Fig. 01 we plot the percentage of the MTA directions 
obtained from the simulated maps that form an angle 
with the observed a dipole direction smaller than a given 
angle (shown on the horizontal axis). The cases of map 
resolutions corresponding to N S id e — 4, 8, 16 are shown. 
The angle between the observed MTA and the observed a 
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FIG. 2: Directions in galactic coordinates for the a (blue) and Dark Energy (green) dipoles, for the Dark Flow direction (red) 
and for the direction of MTA in the 7 years ILC CMB map degraded to N 3 ide = 8 (yellow). The opposite corresponding 
directions are also shown. 




FIG. 3: Percentage of the maximum temperature difference values obtained from the simulated maps bigger than the observed 
maximum temperature difference obtained from the degraded maps with N s id e = 4, 8, 16. 





i n 


6(°) 


MTA (N s id e = 4) 
MTA (N s ide = 8) 
MTA (Nside = 16) 
a dipole 

Dark Energy dipole 
Dark Flow direction 


337.5 ± 14.7 
331.9 ±7.3 
331.9 ±3.7 
320.5 ± 11.8 
309.4 ± 18.0 
282 ± 11 


-9.6 ± 14.7 
-9.6 ±7.3 
-7.2 ±3.7 
-11.7 ± 7.5 
-15.1 ± 11.5 
6±6 



TABLE I: Directions in galactic coordinates for the a [SJ, 
Hfj|and Dark Energy dipoles [Icl|. the Dark Flow and the maxi- 
mum CMB temperature difference (MTA). For the Dark Flow 
direction we have used Ref. Q. The larger scale direction of 
Ref. [f| is consistent with that of Ref. [j but it has signifi- 
cantly larger errorbars. The error on the MTA direction has 
been taken to be equal to the side of the pixel, (12A^ ide ) 
rad. 



dipole direction is indicated by an arrow on each plot. In 
Figs. [SJ [B] we show the corresponding plots where instead 





MTA 


a dipole 


DE dipole 


DF direction 


MTA (N a ide = 8) 
a dipole 
DE dipole 
DF direction 


0.0 
11.4 ± 12 
22.6 ± 18 
52.1 ± 11 


11.4 ± 12 

0.0 
11.3 ± 18 
42.2 ± 11 


22.6 ± 18 

11.3 ± 18 
0.0 

34.4 ± 18 


52.1 ± 11 

42.2 ± 11 
34.4 ± 18 

0.0 



TABLE II: Angular distances in degrees between the Alpha 
and Dark Energy dipoles, the Dark Flow and the MTA di- 
rections. For the MTA direction we have chosen the result 
obtained in the N si d e = 8 case. 



of the a dipole direction we have used the Dark Energy 
dipole and Dark Flow directions respectively. 

The probability that the ACDM simulated maps re- 
produce the observed alignment of cosmic asymmetries 
varies between 1.37% (alignment with a dipole), and 
22.77% (alignment with Dark Flow). The combined 
probability to obtain both a large enough magnitude and 
angular proximity of MTA to the a dipole direction is 
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FIG. 4: Percentage of the maximum temperature difference directions obtained from the simulated maps closer to the a dipole 
than the observed maximum MTA obtained from the degraded maps with Nside = 4, 8, 16. 




10 20 30 40 10 20 30 40 10 20 30 40 

FIG. 5: Percentage of the MTA directions obtained from the simulated maps closer to the Dark Energy dipole than the observed 
MTA directions obtained from the degraded maps with N sil ie = 4, 8, 16. 



shown in Table Mil In particular, the probability to ob- 
tain the observed MTA magnitude (or larger) and the ob- 
served angular proximity to the a dipole direction in the 
context of ACDM varies between 0.5% and 0.1% depend- 
ing on the angular resolution of the WMAP7 ILC map. 
If the probability of obtaining the a dipole magnitude 
in the context of ACDM is also factored in, the proba- 
bility reduces to about one part in 10 7 which is similar 
to the probability for obtaining simultaneously the Dark 
Energy and the a dipoles in the observed directions [l(| ■ 

The last column of Table IIIII is also shown in Table IIVI 
along with the corresponding results for the other two 
anisotropy directions corresponding to the Dark energy 
Dipole and the Dark Flow. 

We stress that the above abnormally low probabili- 
ties assume that the corresponding datasets (Keck+VLT 
quasar absorbers @ , Dark Flow data [H Q , Union2 data 
[39j | and ILC maps [28[) are free of systematic errors. 
The potential validity of these datasets combined with 
the generic nature of the statistical tests applied, assigns 
a particularly low likelihood to the statistical isotropy 
feature of ACDM. 

Nevertheless, the existence of a physical model where 
the alignment of the above axes will appear with a signif- 
icantly larger probability is a prerequisite before putting 
ACDM to disfavor. Even though the qualitative predic- 
tions of Extended Topological Quintessence appear to be 
significantly more consistent with the observed cosmic 
asymmetries than ACDM, a quantitative analysis is re- 
quired before any valid conclusion in favor of Extended 



Topological Quintessence is drawn. Such an analysis is 
currently in progress. 



3. CONCLUSION-OUTLOOK 

We have identified a direction on Maximum Tempera- 
ture Asymmetry (MTA) of the WMAP7 foreground re- 
duced ILC map. Even though the magnitude of this 
asymmetry is consistent with ACDM at the 2a level, 
its direction is abnormally close to other observed cos- 
mic asymmetry axes. The direction of the MTA is close 
to the direction of one of the Cold Spots. This angular 
proximity may imply that this Cold Spot (or the opposite 
located Hot Region) is physically related to the existence 
of other cosmic asymmetry axes. In the context of Ex- 
tended Topological Quintessence, the existence of such 
a feature (Hot or Cold spot) is expected to exist at the 
core of the 'Great Repulser' global defect while in the 
opposite direction an opposite temperature behavior is 
expected. 

The planarity and alignment of the CMB octpulole and 
quadrulole moments may be partly due to a combination 
of two or more features on the preferred plane of these 
moments. Indeed, MTA axis we have identified lies on 
this preferred plane and therefore the MTA may be re- 
lated to the observed quadrupole-octupole alignment [l5[ 

An interesting extension of this project is the deriva- 
tion of the detailed CMB signature predicted by Ex- 
tended Topological Quintessence. Such a derivation 
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FIG. 6: Percentage of the MTA directions obtained from the simulated maps closer to the Dark Flow direction than the 
observed MTA directions obtained from the degraded maps with N 3 ide ~ 4, 8, 16. 





MTA slm s > MTA obs (%) 


@MTA—a ,sims < @MTA-a,obs (%) 


both(%) 


Nside = 4 


16.25 


2.58 


0.48 


Nside = 8 


10.86 


1.37 


0.19 


Nside = 16 


7.12 


1.9 


0.12 



TABLE III: Probabilities of obtaining a simulated CMB map with a maximum temperature difference bigger than the observed 
one and with a MTA direction closer to the a dipole direction than the observed one. 





a(%) 


DE(%) 


DF(%) 


Nside = 4 


0.48 


0.95 


3.37 


Nside = 8 


0.19 


0.52 


2.06 


Nside = 16 


0.12 


0.28 


1.28 



TABLE IV: Probabilities of obtaining a simulated CMB map 
with both a maximum temperature difference bigger than the 
observed one and a MTA direction closer to the a, Dark En- 
ergy dipole and Dark Flow directions than the observed one. 

would involve a cosmological simulation of the evolution 
of the non-minimally coupled 0(3) scalar field that gives 
rise to the recent formation global defect. The linear 
metric perturbations that emerge due to this formation 
can then be numerically calculated [4(3] and the corre- 
sponding ISW effect can be derived in a straightforward 
manner. This numerical analysis can also lead to de- 
tailed predictions about the magnitude and geometry of 
the other cosmic asymmetry axes (Dark Flow, Dark En- 



ergy and a dipole). This analysis is currently in progress. 

Numerical Analysis Files: The data, Mathematica 
and Healpix program files used for the numerical anal- 
ysis files may be downloaded from |http : //leandrosT] 
[physics .uoi .gr/mta| 
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